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SYMBOLS 


State (Ns,Ns) or Output (No,No) Weighting Matrix 
Control (Ne, Nc) Weaghting Matrix 

Contre Gain Metrix. (Nc, Ns) 

Control (No,Nc) or Noise (No,Ng) Feedforward 
Matrix 

Open-Loop Dynamics Matrix (Ns,Ns) 

Control Distribution Matrix (Ns, Nc) 

State ‘D¥sSturbance Distribution Matrix’ -(Ns,Ng) 
Measurement Scaling Matrix (No,Ns) 

Estimator Gain Matrix (Ns, Nob) 

Number of Ccntrols 

Number of Process Noise Sources 

Number of States 

Number of Okservations or Measurements 

White Process Noise Covariance Matrix (Ng,Ng) 
White Meas. Noise Covariance Matrix (No,No) 
Steady-State Covariance Matrix of Control (Nc,Nc) 
ContEcl! Vector! (Nc ,'1) 

White Measurement Noise Vector (No,1), with 
Zero Mean and Covariance Matrix R 

White Process Noise Vector (Ng,1), with 

Zero Mean and Covariance Matrix Q 

Constant Disturbance Vector (Ng, 1) 

State Vector (Ns,1) 

Derivative of State Vector (Ns,1) 

Estimate of State’ Vector’ (Ns, 1) 

Derivative of Estimate of State Vector (Ns,1) 
Output Vec ter’ (No 74) 


Measurement Vector (No, 1) 


I. INTRODUCTION 

The purpose,,of this) thesis 4 is) to (deserzbe, fie medi ited] 
tion and additions: te.the, existing, FORTRAN, program g(OPTSyo) 
which is used in the study and) application of ») Oftinal 
Systems Ccntrol theory. 

The OPtimal SYstems control program was originally 
developed by Halil -[Ref.. U]}) to. Support» Diss research jan 
rotary-wing aircraft control systems. Later program modifi- 
cations were made by Walker [Ref. 2] and Liu [Ref. 3] of 
Stanford University, and are designated OPTSYS 4 and OPTSYS 
5 respectively. OPTSYS modifications made by Hoden [Ref. 4] 
were primarily devoted to creating a user-friendly interac- 
tive version (OPTSYSX) of the OPTSYS 4 progran. The latest 
modifications by Diel [Ref. 5] allowed the user to save 
matrices for subsequent runs of the OPTSYSX program, and 
formed data file sets for the time response program OPTCALC. 

The intent of this thesis work was to develop an inter- 
active program to plot the Pole-Zero map, Root-Locus, Bode, 
Nyquist, and Nichols responses to the open loop, closed loop 
noise, and; compensator transfer) ~functsions-4it wa State 
Variable Control System which has been developed using 
OPTSYS Program. Mincr modifications to the OPTSYSX Program 
were necessary to ccnstruct, aj sdataifilessetsip for plotting 
the Pole-Zero Map, andin calculations for the Root-Locus, 
Bode, Nyguist, and Nichols plots. 

It is assumed that the user is familiar with the kasic 
concepts cf Control Theory and Optimal Systems Design. The 
symbol/naming conventions of Bryson [Ref. 6]. are used in 
the discussion of program operations and descriptions of 
problems using the OFTISYS System. 


An overview of the OPTSYSX Program capabilities and of 
modtiacations tos thevexisting programs is presented first. 
This is followed by a description of the program (OPTGRAPH), 
which was developed tc plot the Pole-Zero Map and to perform 
calculations for the Root-Locus, Bode, Nyquist, ‘and Nichols 
pLots. 

This work concludes with examples of various types of 
problems demonstrated in the interactive mode, including a 
copy of each terminal session with the final results. 
Complete program  astings,, for the OPTSYS EXEC. program, 
CPTSYSX progran, and CPTGRAPH program are included in 


appendices A,B, and C respectively. 
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II. THE OPTSYSX COMPUTER PROGRAM 


A. BACKGROUND 


OFTSYSX is a  double-precision, interactive FORTRAN 
program employing modern control theory analysis techniques. 
Its “capabilities “include the calculation of the cpen-icop 
eigensystem, and the stationary closed loop system; the 
synthesis of regulators and filters; along with pcwer spec— 
tral density, and modal distribution computations. The 
modifications introduced to the’ OPISYSX "program by “Ehts 
thesis work do nct affect the program's original 


capabilities. 


Be OVERVIEW 


OFTSYSX is an extremely large and complex program with 
over 3000 lines of ccde. To use the program dn. tts small 
version (dimensioned for .a 32. X% 32, "%E", “gr, and =~"H" 
Matrices), the user must extend his Virtual Machine (VM) 
memory capacity beyond the default VM memory capacity cf 720 
kilobytes to 1024 kilobytes, 

Any significant increase in the OPTSYSX program size, 
with the the resulting requirement for additional memory 
capacity, would cause user difficulties. For this reason 
the task of obtaining the Pole-Zero maps, Root-Locus, and 
Bode, Nyquist, Nichols plots for the open loop, noise and 
compensator transfer functions was relegated to a separate 
program (OPTGRAPH). ‘Three data files (OPTGROL DATA, OPTGRNO 
DATA, and OPTGRCM DATA), containing the open loop, noise, 
and ccmpensator transfer functions, are formed by the OPTSYS 
program to be transfered to the OPTGRAPH program. The 


OPTSYS EXEC program provides an interface between the 


We 


OPrSy SX ‘progran; the time response programs (OPTCALC and 
OPTPLCT), and the OPTGRAPH progran. 


C.2) OPTS YSX®* MODIFICATIONS 


QETSYSX program modifications consist of the addition of 
write statements, three flags, and a short routine to inter- 
pect user inputs. 

Write statements,to input data to the OPTGROL, OPTGRNO, 
and CPTGRCM DATA files were added to the main progran, 
subroutine INNER, subroutine ZEROS, and subroutine POLES. 
System information ccnsisting of number of states (Ns), 
number of controls (Nc), number of measurements (No), number 
of process noise scurces (Ng), type transfer function 
(ITFX), Markov Parameter (IE), and two flags to be used by 
the CFTSYS EXEC and the OPTGRAPH programs is inputed from 
the main program. The write statements added to the subrou- 
tine INNER, and Subroutine POLES input the poles for the 
open locp and noise transfer functions, and compensator 
transfer function respectively. The zeros, the numerator 
order, the gain, and the input and output numbers for all 
three transfer functions are obtained from the subroutine 
ZEROS. 

The first of the two flags, input with the system infcr- 
mation from the main frogram to the three data files, serves 
as dap Marker Lor the CPTSYS/EXEC*°to locate’ the-correct data 
line for reading system inputs. The second flag, which is 
set by the routine tc interpret user inputs, serves asa 
Sigidia: fOr the OPTSYS SEXEC3tO! either!’ load the® OPTGRAPH 
pPEeogram if.any one) o£+the* three? transfer functions” was 
calcudated’ by’) the OFISYSX-progran, and abort loading the 
OPTGRAPH program if none of the three transfer functions 
were calculated. ThePthirdsutlaga (il rx) is an existing 
OPTSYSX frogram flag. The flag (ITFX) is passed from the 
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subroutine TF to the subroutine ZEROS to identify the calcu- 
lation of either an open sloepygsmorse, or compensator 
transfer function zeros. 

The short routine added to the main program of OPTSYSX 
reads the’ flags ITF, 2TE2, I9F3,, and RET. The flag sent 
to the OPTGROL, OPTGRNO, or OPTGRCM DATA files is set to 
allow loading of transfer function data file to the OPTGRAPH 
program if the user selects either of the transfer function 
options number two (calculate poles, residues, and zeros - 
TRE ="), 29R2Z. =e Obes veres.— ei) humber three (calculate 
poles, and zeros = ITFI = 2, .1ITF2 =2, or, iTF3)=)2) wea 
calculating the open loop, noise, or compensator transfer 
FUNCL TONS. The selection of any other transfer function 
option will not calculate the transfer function zeros. The 
flag IRET is read to determine if the user desires to rerun 
a problem without exiting the OPTSYSX progran. Should the 
user elect to rerun a problem without exiting OPTSYSX (IRET 
= 1), the three data files are set back to the beginning 
(REWIND) to accept the new problem system data. 


De OFTSYS EXEC MODIFICATIONS 


The OPTSYS EXEC was written by Diel [Rek. 5] tc, anter— 
face between the OPISYSX program andthe time response 
programs OPTCALC and CPTPLOT. The modifications to the exec 
program were to .add the option:)to ~select running, the 
OPTGRAPH program to the exec menues and a_ short section to 
interface the OPTGRAPH program with the OPTSYSX program. 

The added section to interface between the OPTSYSX and 
OPTGRAPH programs defined loader size (LDRTBLE), file defi- 
nitions.- (FILEDEF);,, and, ,;textslibraries s{TKTLE Bjgttors! the 
CPTGRAPH program. The first line of the three data files is 
read to ensure that there has been at least one of the 
three transfer functions calculated by the OPTSYSX prcgram 
prior to loading OPTGRAPH progran. 


a3 


THe OPISYS MEKECVrOULIne@isjaliso used to send plots to 
tHCee VERSATEC? Splottecr “trom thes nOPTPLOT program and’ the 
OPTGRAPH progran. 
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III. THE OPGRAPH PROGRAM 


A. PROGRAM OVERVIEW 


OPTGRAPH iS an interactive FORTRAN program which uses 
the transfer function poles and zeros calculated by OPTSYSX 
to calculate, tabulate, and plot the Pole-Zereo Map,. Eode;, 
Nyguist, Nichols, and Root-Locus responses for either the 
open loop, closed loop noise, or compensator transfer 


fuNnceElons. 
1. Frogram Language 


CFTGRAPH is programmed in FORTRAN following the 
conventions of the IBM System /360/370/ FORTRAN IV language. 
CPTGRAPH has been ccmpiled and run under both FORTRAN IV 
(G1) and FORTRAN H (extended) compilers on the IBM 3033. 


2. Graphics Package 


This program uses the Display Integrated Software 
System and Plotting System Software (DISSPLA) developed and 
distributed by the Integrated Software Systems Corporation 
(ISSCC) of “San. Diego, Calviornma. The DISSPLA package is a 
library cf FORTRAN subroutines using FORTRAN IV conventions. 


3. ILikrary Subroutines 


CPTGRAPH uses the International Mathematical and 
Statistical Library (IMSL) ysubroutimesZRPOLY to recalculate 
the roots of a polynomial equation. ZRPOLY iS. a deubike 
precision FORTRAN subroutine to which is input a polynomial 
with real coefficients (double precision) ine terms xo 
decreasing powers of the polynomial. The subroutine finds 


the rects of the polynomial and returns a double precision 


cS 


complex array containing the roots. GRPOLY is Sicapab le ot 
accepting polynomials greater than 0 order and less than 101 
order. 


CPTGRAPH has one main program and 24 subroutines. 
The main program and its subroutines may be divided into 
four kasic categories: 
T) -"File’ Data Input 
2) Interactive Data Input 
3) Plot Setup and Sequencing 
Ly. -Calculation 
A brief and general description of the program and 
its subroutines will follow in subseguent sections. 


Be. GENERAL PROGRAM OPERATION 


OFTGRAPH was written to satisfy two specific objectives. 
Pic ftest 1S £O Create a program to be used as a instruc- 
tional tools £or-students» taking ‘controls -rcelated, courses, 
and the second is as method to assist in the analysis) Of 
actual large order ccntrol system problems. 

To accomplish these objectives the emphasis while 
wEEting, this) {program “was ito Make; the, program asi) user 
oriented as possible Fy attempting to eliminate ambiguities 
and providing features which would protect the user from 
inadvertent wrong entries. An effort was also made to mini- 
mize the Virtual Machine (VM) memory required to execute the 


program to facilitate the analysis of a large order systen. 


Le User Protecticn Features 


Three basic methods were used to protect the user 
from an abnormal program termination in the case of an inad- 


WErtent wrong entry. Tic Eotewass tO) construct, the 
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subroutines to display a summery of the user inputs and 
llowing changes prior to leaving the subroutine. For the 
case where the user input is expected to fall within a given 
range, the input is examined and the user is either allowed 
to continue for an “input within the range or issited an) eeron 
or warning and reasked the question for an input outside the 
range. The final method is used by the three subroutines 
which read the user input from the screen. The subroutines 
RDINT, RDREAL,and RICHAR expect an integer input, real 
humber input, and a lcgical "yes" or "no" respectively. In 
the event that a "null" line is inadvertently entered, the 
user is issued a warning and allowed another opportunity to 
enter the the correct input before abnormally terminating 
the program. The entry,.of .two "null" lines also offers =the 
user the option to to exit the OPTGRAPH program at cther 
than normal program exit points. The subroutine RDCHAR also 
examines the input and issues a warning if the input is not 


e1ther “yes or “not. 
2. iarge Order System 


The large order control problem being considered for 
analysis by the OPTSYSX and the OPTGRAPH programs is the 
X-29A aircraft longitudinal axis backup mode systen. This 
system has a “(98" x 968) °F" matrix , “a (2 x 99) TH! Me thal ee 
and a S(98 X 1) "CG matrix. 

To accommodate this large order system the OPTGRAPH 
program was dimensioned to accept a system with a maximum of 
99 states “(Ns), ‘12 ‘controls’ (Nc), 2 process noise sources 
(Ng), and 12 measurements (No). The Virtual Machine (VM) 
memory requirements for the program were minimized by the 
reuse of memory Locations allocated “for Storage. Of Calci- 
lated data after that data has been either tabulated or 
plotted. This methcd of memory reuse has allowed the 
CPTGRAPH program to ke dimensioned to accommodate a large 
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order system and still operate with less than the 1024 
RETORYIES Vi Stmenosy: Tequirede tolrin? * OPTSYSX*®in its small 


version. 


Cw. SYSTEM/MODEL DESCKIPTION 


The system equations used by the OPTSYSX program for the 
calculation of the transfer functions are of the state vari- 
able form. The system equations are: 
system model 

xdot = [F]*x + [G]*u + [GAM )]*w Gin 
measurement equation 

Z=— PH j*x + (Dieu + y¥ (Bez) 
estimator equation 

LedGte— ser \oce: fC jaa ite Oi | *i(z es ee] *xie)s 6 £33) 
open loop transfer function 

Lage )s( ap (slate (3.4) 
closed lcop noise transfer function 

Ph hs (Be SeE hl +44 Peany (3-5) 


compensator transfer function from measurement to input 


beI* (st RISE eLG j* caer xe aaie ty Sek | (3.6) 


where 


u = Control vector (Ns X ‘1) 
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W = White process noise vector (Ng X 1) 


xX = State vector (Ns X 1) 

xdot = Derivative of the state vector (NS°x 1) 
y = Output ®vector (No x71) 

Zz = measurement vector (No X 1) 


{C] = control gain matrix (Nc,Ns) 

{D] = control feed-forward distribution matrix (Nc, Nc) 
of the controlevector? (u) 

[F] = open-loor dynamics matrix (Ns, Ns) 
(System matrix or Plant matrix) 

[G] = controlydustributioninatrmx (Ns, Ne) 

[GAM] = state disturbance distribution matrix (Ns,Ngqg) 
{[H1] = measurement distribution matrix (No, Ns) 


[I] = identity matrix (Ns,Ns) 


D. PROGRAM ORGANIZATION 
1. Main Program 


The main program presents a menu allowing the user 
to select for analysis either the open loop transfer func- 
tion, the closed loop noise transfer function, the compen- 
sator transfer functicn, or exit the OPTGRAPH program. The 
data file for the selected transfer function is examined by 
the main program to ensure that the file contains the 
required transfer function data. Should the data file be 
incomplete the program gives the user a warning message 
indicating the nature of the missing data and allows the 
user tc select another option. 


The main program menu: 
OFIGRAPH 


DO YOU DESIRE TO ANALYZE: 
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1. OPEN LOOF TRANSFER FUNCTION 
2.) NOISE (ERANSFER FUNCTION 
3. COMPENSATOR TRANSFER FUNCTION 
4. EXIT OPTGRAPH 

ENTER OPTION NUMBER. 


Data is read from the selected file by the main 
program and formed into three arrays and seven integer vari- 
ables. The exchange of data between the main program and 
the subroutines, and between the subroutines is by the 
extensive use of labeled common statements. The data trans- 
fered through the subroutine call statements is limited 
Permmarily to Elags, and constants. The program moves to its 
primary sequencing subroutine (GRAPH) upon the completicn of 
file data entry. 


2. FErogram Sequencing 


The subroutine GRAPH first presents the user witha 
Summary cf the the transfer system received from the data 
file. The user is then asked which transfer function he 
desires to to analyze using the conventions of the OPTSYSX 
program describing the transfer functions by input number 
and output number. A menu is presented which allows the 


Msen, to select graphical dnd/or tabular system response 


(Pole-Zero map, Root-locus, Bode, Nyquist, or Nachols) | 
select another transfer function, Or exit the OPTGRAPH 
progran. Upon completion (of a graphical and/or takular 


response option, or change of transfer function option, the 
peegram LTeturns to the GRAPH menu to allow the user to 
choose another option or exit to the main program. 


Subroutine GRAPH data summary and menu: 
OFTGRAPH 
THE OPEN LOOP TRANSFER SYSTEM OBTAINED FROM OPTSYS CONTAINS: 


i, SUSBER OF STATES = 4 
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2. NUMBER OF CCNTROLS (INPUTS) = 1 
3. NUMBER OF MEASUREMENTS (OUTPUTS) “= 
4. MARKOV PARAMETER = 10.0**- 6 


CLEAR SCREEN TO CONTINUE 


OPTGRAFH 


AN OPEN LOOP POLE-ZEKO, ROOT LOCUS, BODE, NYQUIST, AND/ 
OR NICHOLS PLOT IS DESIRED FOR: 
INEUL: #2 


CUTPUT #2 


OPTGRAFH 


OPEN LOOP TRANSFER FUNCTION 
INPUT # = 1 
OQUTEUT # 1 
DO YOU DESIRE GRAPHICAL RESPONSE AND/OR TABULAR DATA FOR: 
1. POLE-ZERO MAF 
2. ROOT-LOCUS 
3. BODE 
4. NYQUIST (POLAR PLOT) 
5. NICHOLS 
6. ANOTHER INPUT/OUTPUT COMBINATION (SAME TRANSFER 
FUNCTION) 
7. ANOTHER TYPE TRANSFER FUNCTION /EXIT OPTGRAPH 
ENTER OPTICN NUMBER. 


3. Graphic and/or Tabular System Response 


The selection of a graphical and/or tabular system 
response option for an Pole-Zero map, Root-Locus, Bode, open 
loop Nyquist, or ofen loop Nichols calls the subroutines 
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PEERO RW ROLO,MBODE, NYCSDTZS0r, NECHOLO respectively.) While the 
response of these subroutines is different, their crganiza- 
ton isi:suma lars 

The Subroutines! PZERO,YRTILO,YS BODE; ONYOST,- and® NICHOL 
primarily act as sequencing and plot setup subroutines. The 
nrerachive sy questioringumof ithe wser ©for fgraphic/tabular 
response information and the system response calculation has 
been relegated to other subroutines. The subroutines PZERO, 
wet tOrny SOME, ~NYOST weandsNICHONG are dividedfinte twot major 
sections with the first being devoted to producing a graph- 
ical output and the second to producing a tabular output. 


The user is presented a subroutine menu offering the 


cption fcr graphic response, tabular response, or exiting 
the subroutine. The program moves to either the first 
section (graphic response), second section (tabular 


response) or back to the subroutine GRAPH menu depending on 
the option selection. The program returns to the subroutine 
menu upon completion cf the of either the graphic response, 
or takular response section, allowing the user to select 
another response opticn, or exit the subroutine. 

Features that are common to the five subroutines 
will ke described here and features that are unique to one 
of the subroutines will be described in the follcwing 
section under the apprepriate response heading. 

For a graphic response, the user is given tke choice 
G ewe Printer opticens (TEK 618, . or. VERSATEC). With the 
selection of the TEK 618, the plot page size is defaulted to 
lm@anenc= by 8.5 imches, and the plot will be presented on 
the TEK 618 screen. The VERSATEC option gives the user the 
OpELon te create up to a 21] inch by 21 ianch. plot. (maximum 
EGE Ene VERSALEC printer) and causes a DISSPLA METAFILE to 
be created. A scalang subroutine (PSCALE) scales the plot 
heading, legend, etc. im proportion to the page width 
selected. The VERSATEC option offers the the advantage of a 
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high quality print; fromy(the* 7VERSATEC printer) usingdtue 
VERSATEC printer, option an’ the, OPTSYS) ExECse > Thelereadtvong of 
a DISSPLA METAFILE “also “gives the wser |i(efterm he yexases 
optsys) the opportunity to use printers other than the TEK 
618 or->o VERSATEC. A major disadvantage with the VERSATEC 
ception is that the user must exit the OPTGRAPH program 
kefore he can print a graphic system response. For) this 
reason, it is recommended that the system response first be 
plotted: on thes! TERK, 618 prior itouaselecting thea VERSAIZEe 
option. 

For a tabular response,the user is presented a menu 
which gives him the choice of three devices (screen, 
printer, or disk) to send the tabular output. The user is 
asked if he desires to make any changes, after the tabular 
data has been calculated. If the answer is affirmative, he 
is presented a menu and allowed to make changes and rerun 
the tabular data routine. 

Examples of these features are included in the 
following sections with the interactive examples for the 


system responses. 


E. INTERACTIVE EXAHFIES 


The open-loop, closed loop noise, and compensator 
transfer functions are are handled identically by the 
OPTGRAPH program, with the only difference being the program 
and graph headings identifying the transfer function tyre. 
The frequency response programs for the Bode, Nyquist, and 
Nichols responses use identical setup and tabular data 
routines. 

A gocd overall representation of the OPTGRAPH program 
capabilities will be demonstrated by single interactive 
terminal session examples for the pole-zero, root locus, and 


bode resfonses using selected transfer functions. Tabular 
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data examples will ke provided for the Pole-Zero, Root 
Locus, and Bode responses, and graphic response examples 
will be provided the Fole-Zero, Root Locus, Bode, Nyquist, 
and Nichcls responses. 


F. CHARACTERISTIC EQUATION ANALYSIS 


Phe subroutines PZERO and RTLO are similar in their 
basic formats. The pole and zero locations for the selected 
EbanSter Lunction input and output numbers are read into two 
Single dimension arrays. The extraneous zeros, calculated 
by the OPTSYSX progtran, are eliminated by comparing the 
zeros to the Markov parameter sent with the system infcrma- 
tion from the OPTSYSX progran. For a graphical response, 
the user is asked to define the plot limits in terms of 
X-coordinates (real axis) and Y-coordinates (imaginary 
axiS) « Points for the plot which fall outside these limits 
will be ignored. System information (tranfer function 
input number, output number, and (DC) gain) is listed in the 


graph legend. 
1. Pole-Zero Map 


The subroutine PZERO converts the double precision 
humbers for pole and zero locations to single precision for 
compatibility with the with the Graphics Package and flots 
the poles and zeros within the plot limits defined by the 
user for the graphical response routine. For tabular data 
response the user is given system information (numerator 
order, denominator order, and transfer function (DC) gain) 
for the selected transfer function input number and outrut 


humber, and the pole and zero locations. 
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2. Fole-Zero Map (Interactive Example) 


The following example of a simplified open loop transfer 
function for an “airplane with an autopilot, 2n the Jeng ru. 
dinal mode was obtained from [Ref. 7,p.315]. 


The aircraft system is represented by: 


ee a 


Onatle, oSep 0. of 
Org iOfetye alicn “Os. wil 
| 
| 


CF] 
Qe 940s Oey elie 
Oi. 16..—= hae -3. 


CG] =f iG- 


— 
—_ 
e 

— — — — 


[GAM] = [0] 


fH} =f “los Oe Oem OG. 


a. Example cf Graphic Response (Terminal Session) 


BEGIN RECORDING OF TERMINAL SESSION 
Rs T=O..09 7050201 S207 
EXECUTION BEGINS... 


OPTGRAPH 
DO YOU DESIRE TO ANALYZE: 
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1. OPEN LOOP TRANSFER FUNCTION 
2. NOISE TRANSFER FUNCTION 
3. COMPENSATOR TRANSFER FUNCTION 
4. EXIT OPTGRAPH 
ENTER OPTION NUMBER. | 
4 


OPTGRAPH 


LHE, OPEN TOCERIRANSEPER) SYSTEM OBTAINED FROM OPTSYS 
CONTAINS: 
1.. NUMBER OF STATES = 4 
Z2--) NUMBER OF SCONTROLS: (TE NPUTS) = 1 
3. NGMBER OF MEASUREMENTS (OUTPUTS) 
i, UARKOV PARAMETER = 10.0**— 6 


i} 
— 


CLEAR SCREEN TO CONTINUE 


OPTGRAPH 


KN OPEN LOOP POLE-ZERO, ROOT LOCUS, BODE, NYQUIST, 
AND? OR? NICHOLS. PLOT=~ IS DESTRED” FOR: 
ENPUY 2-2 


OUTPUT #2 


OPTGRAPH 


OPEN LOOF TRANSFER FUNCTION 

INPUL # = 1 

OUTEUT # 1 
DO YOU DESIRE GRAPHICAL RESPONSE AND/OR TAEULAR 
DATA FOR: 

12) FOLE-ZERO” MAP 
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2.  ROOT-LOCUS 
3. EODE 
4. NYQUIST (POLAR PLOT) 
5. KICHOLS 
6. ANOTHER INPUT/OUTPUT COMBINATION (SAME 
TRANSFER FUNCTION) 
7. ANOTHER TYPE TRANSFER FUNCTION /EXIT 
CPTGRAPH 
ENTER OPTICN NUMBER. 
: 


POLE=Z ERGs MAE 

OPEN LOOP TRANSFER FUNCTION 
DO YOU DESIRE: 

1. GRAPHICAL OUTPUT 

2. TABULAR DATA 

35. GUIT SUBROUTINE 
ENTER OPTICN NUMBER. 
1 


POLE-ZERO MAP 
OPEN LOOP TRANSFER FUNCTION 
PLOTTER SELECTION ANDSP AGE SIZE 
(NOTE: PAGE TS 117.0 ZNEHES BY 825 INCHES Wii 
SELECTICN OF DEKG 18) 
WHICH PLOTTER DO YOU DESIRE: 
1. TEK618 
2. VERSATEC 
CHOGSE OPTION 19 0OR 2 


PAGE SIZE (MAXIMUM = 212.0 INCHES BY ©2 1-00 


INCHES) 
HEIGHT = 


Pail 


WIDTH, = 
8 
PLOTTING LIMITS’ PORY GRAPH 
X AXIS (REAL AXIS) 
xX AINEIMUN = ? 


=) 5 
X MAXIMUM = ? 

5 

Y AXIS (IMAGINARY AXIS) 

Yo PEN = 2 

5) 
Y MAXIMUM = ? 

5 


PLOTTING IIMITS FOR GRAPH 
Xe AxaS (REAL AXTS) 


X MINIMUM = = Ss O10 
X MAXIMUM = 5.00 
Y AXIS (IMAGINARY AXIS) 
Y MINIMUM = ='54,0'10 
Y MAXIMUM = 5-00 


DO YOU“ DESTRE TO MAKE ANY -CHANGES ? 
TYP Reese OR NOt. 


Nn 
PLOT HEADING 


HOW MANY LINES OF HEADING DO YOU DESIRE ? 
( 3MAX ) 
ENTER NUMBER OF LINES. 
(ENTER "0" FOR NO HEADING) 


Ay MAXEMUN" OF 32 CHARACTERS: PER-LINE IS 
ALLCWED 
LENE 1S 

auneratt with an autopilot in 
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LINES2 ES. 

the longitudinal mode 
LINE: Ses: 

ref. 73 pagerats 


PLOT HEALING 


LINE. (2 ATRCRAFT WITH AN AUTOPILOT =F 
TNE os THE LONGITUDINAL MODE 
LINE, 35's REF. 73 PAGE 315 


DO YOU DESIRE TO MAKE ANY CHANGES 2 
TYPE YES" OR UNO. 
nN 
>USING A PRE-ALLOCATED DATASET FOR UNIT FT17F001. 
>USING A FRE-ALLOCATED DATASET FOR UNIT FIIGFOOR 


A DISSPLA METAFILE HAS BEEN CREATED 
CLEAR SCREEN TO CONTINUE 


POLE-ZERO MAP 
OPEN LOOP TRANSFER FUNCTION 
DO YOU DESIRE TO MAKE ANY CHANGES TO: 
1. PLOTTER “7 PACE PSEZE 
2- GRAPH Eauirs 
3. HEADING 
4. NO ADDITIONAL CHANGES - PLOT POLE-ZERO 


MAP 
5. NO CHANGES —- EXIT POLE-ZERO PLOTTING 
ROUTINE 
ENTER OPTION NUMBER. 


5 


POLE~ZERO MAP 
OPEN LOOP TRANSFER FUNCTION 
DO “YOU DESIRE 
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1. GRAPHICAL OUTPUT 
2- TABULAR DATA 
S-2SOULTISUBROUTINE 
ENTER OPTICN NUMBER. 
3 


OPTGRAPH 
OPEN LOOP TRANSFER FUNCTION 


INPUT @iunse al 

OUTPUT» #5= 1 
DO YOU DESIRE GRAPHICAL RESPONSE AND/OR TABULAR 
DATA FOR: 

1. FOLE-ZERO MAP 

25, — HOOT-LOCUS 

Ze. BODE 

4. NYQUIST (POLAR PLOT) 

5- NICHOLS 

6. ANOTHER INPUT/OUTPUT COMBINATION (SAME 

TRANSFER FUNCTION) 
da) ANOTHER YRS, TRANSFER, FUNCTION EXIT 


CPTGRAPH 
ENTER OPTICN NUMBER. 
7 
OPTGRAPH 


DO! YOU) DIPSTRET TOM ANALYZE: 
To) CLEN LOOP TRANSFER FUNCTION 
2- NOISE TRANSFER FUNCTION 
3. COMPENSATOR TRANSFER FUNCTION 
4. EXIT OPTGRAPH 

ENTER OPTION NUMBER. 

4 
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Rs T=2. 8074-6" 4492 456415 
record off 
END RECORDING OF TERMINAL SESSION 


The preceding example follows at the end of the 


chapter as figure: sen. 
b. Example of Tabular Data (Terminal Session) 


The entering and exiting option menus for the 
tabular data terminal session are identical to the graphic 
response terminal session option menus, and have been elimi- 
hated frcem the follcwing terminal session. 

BEGIN RECORDING OF TERMINAL SESSION 
R3; T=0.0170.02) 19351203 
EXECUTION BEGINS... 


POLE=ZERO MAE 

OPEN LCOCF TRANSFER FUNCTION 
DO YOU DESIRE. 

i. GRAPHICAL OUTPUT 

2- TABULAR DATA 

36 QUIT SUBROUTINE 
ENTER CPTIICN NUMBER. 
2 


POLE-ZERO MAP 
OPEN LOCP TRANSFER FUNCTION 
DO YOU DESIRE TABULAR OUTPUL=1OF GO. TO? 
1. SCREEN 
Zs PRINTER 
3. DISK (OPGRAPHSLESTING) 
ENTER CPTICN NUMBER. 
1 
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POLE-ZERO MAP 
OPEN LOOP TRANSFER FUNCTION 


INEUT NUMBER = 1 
CUTFUT NUMBER = 1 
TRANSFER FUNCTIEN (DC) GAIN = 0..1000D+01 
DENCMINATOR ORDER = ~ 
POLE LOCATIONS 
REAL PART IMAGINARY PART 
0.0000D+00 0.0000D+00 
0.1000D+01 0.0000D+00 
=0520,00D+ 01 0.3464D+01 
-0.2000D+01 -0.3464D+01 
NUMERATOR ORDER = 1 
ZERO LOCATIONS 
REAL PART IMAGINARY PART 
-0.1000D+01 0.0000D+00 


DO FOU, DESTRE, TO, CHANGE) OUTPUT, DEVICE? 
eer BbYES. “SOR, o* NOM. 


POLE-ZERO MAP 
OPEN LOOP TRANSFER FUNCTION 


DO YOU DESIRE: 


1. GRAPHICAL OUTPUT 


2- ##TABULAR DATA 
34 OULT SUBROUTINE 


ENTER CPTIICN NUMBER. 
3 


Bee t= 1460/3.23 2953300 


record off 
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END RECORDING OF TERMINAL SESSION 


3. Roct-Locus 


The subroutine RTLO creates a closed loop system 
(with negative unity feed-back) from the selected transfer 
function input numbers and output numbers. Transrer fune— 
tion poles and zeros are formed into two separate pclyno- 
Mials by the subroutine MAKPOL. The numerator polynomial 
(zeros) is multiplied by the gains, and the two polynomials 
are added to form a single polynomial. Complex roots of the 
of the pclynomial are calculated? by the IMS’ Yebracy sukrcou- 
tine ZRPCLY. The ccmplex roots are separated into real and 
imaginary components for either plotting or tabulating. 

The plotting routine first plots the transfer tunes 
tion pole and zero lecations then plots the root tecatren— 
as they are received from the subroutine ZRPOLY. For the 
plotting routine the gain interval specified by the user is 
divided into 2000 evenly spaced points for the subroutine 
ZRPOLY, to caleulate saces. The poles, zeros, and roots are 
plotted within the plct limits specified by the user. 

For the tabular data routine the user is given the 
ception tc choose between 1 and 500 points for the sukroutine 
ZRPOLY to calculate roots. The tabular data output provided 
for the user is system information (numerator order, denoni- 
nator order, and transfer function (DC) gain)--“for’* thie 
selected transfer function input number and output number, 
the pele and zero locations, and the root locus gains and 
roots 

A positive unity feedback system may be analyzed by 


specifying a negative gain range vice a positive gain range. 
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4. Root-Locus (Interactive Example) 


The following interactive example is of the ofen 
Poop. tEanster Lunctien for an aircraft with an autepilot in 
the longitudinal mode described in Pole-Zero interactive 


secticn. 
a. Example of Graphic Response (Terminal Session) 


BEGEN RECORDING OF TERMINAL SESSION 
Roel —0. 0170.02 20228337 
FLECUTION BEGINS... 


OPTGRAPH 


EOS TCU DESTRE TO ANALYZES 
i. OPEN LOOP TRANSFER FUNCTION 
Za NOUSE TRANSFER FUNCTION 
aS COMPENSATOR TRANSFER FUNCTION 
oe Exit ORTGRAEH 

ENTER CPILICN NUMBER. 

1 


OPTGRAFH 


Pee CreEN LOOP TRANSFER SYSTEM OBTAINED FROM OPTSYS 


CONTAINS: 
1. NUMBER OF STATES = 4 
Zs NUMBER OF CCNTROLS (2NPUTS) <= 
3. NUMBER OF MEASUREMENTS (OUTPUTS) = 1 
4. MARKOV PARAMETER = 10.0**- 6 


CLEAR SCREEN TO CONTINUE 


OPTGRAFH 
AN OPEN LOOP POLE-ZERO, ROOT LOCUS, BODE, NYQUIST, 
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AND/ORK NICHOLS PLOT 25 DESTREDVFOR: 
INPUT # ? 


OQULRUT # @ 


OPTGRAFH 


OPEN LOCE TRANSFER FUNCTION 

IWPUT FOS] 4 

CUTEUT # 1 
DO YOU DESIRE GRAPHICAL RESPONSE AND/OR TABULAR 
DATA FOR: 

1.  POLE-ZERO MAP 

2.  FOCT-LOCUS 

3. BODE 

4. NYQUIST (POLAR PLOT) 

5. NICHOLS 

6. ANOTHER INPUT/OUTPUT COMBINATION (SAME 

TRANSFER FUNCTION) 
7. ANOTHER TYPE TRANSFER FUNCTION /EXIT 


OPTGRAPH 
ENTER CPTION NUMBER. 
2 

ROOT-LOCUS 


OPEN LOCE TRANSFER FUNCTION 
LO) YOU DESZTRES 
1. GRAPHICAL OUTPUT 
2. TABULAR DATA 
3. QUIT SUBROUTINE 
ENTER CPTION NUMBER. 
1 
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ROOCT=LOCUS 
OPEN LOOE TRANSFER FUNCTION 
EL@TLER SELECTION ANE PAGE SIZE 
eles PACE JS 1150 -2NCHES BY -c.5 ENCHES WITH 
SELECTICGN OF TEKG 18) 


WHICH PLOTTER DO YOU DESIRE: 
1. TEK618 
2. VERSATEC 

CHCCSE OPTION 1 CR 2 


2 
PAGE SIZE (MAXIMUM = 21.0 INCHES BY 21.0 
INCHES) 
HEIGHT = 
6 
WIDTH = 
8 


PEOTTING LIMITS FOR CRAPH 
X AXIS (REAL AXIS) 
X MINIAUMC =" ? 


aS 
X MAXIMUM = ? 
5 
Y AXIS (IMAGINARY AXIS) 
Y MINIMUM = ? 
=, 
VY oMAXEMUN c= 0? 
5 


PLOTTING LIMITS FOR GRAPH 
X AXIS (REAL AXIS) 
X MINIMUM -15.00 
X MAXIMUM 5.00 
Y AXIS (IMAGINARY AXIS) 


iT} 
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Y ALVES -5.00 
Y MAXIMUM 5.00 
DO YCU DESIRE TO MAKE ANY CHANGES = 
TYPE VYESM (OR) NOW. 


GAIN KANGE FOR ROOT=LECUS PLOT 
MINIMUM GAIN & 


tl 
Pe) 


MAXIMUM GAIN 
500 


GAIN RANGE FOR ROOT-LCCUS PLOT 
MINIMUM GAIN 0.000D+00 
MAXIMUM GAIN = 0.500D+03 
DO YOU DESIRE TO MAKE ANY CHANGES ? 
TYPE CYES” OR VUNG: 


ELOT HEALING 
HOW MANY LINES OF HEADING DO YOU DESIRE ? 
( 3 MAX ) 
ENTER NUMBER OF LINES. 
(ENTER "O" FOR NO HEADING) 


A MAXIMUM OF 32 CHARACTERS PER LINE 5S 
ALLOWED 
LINE OU CLS 
aircraft wath an autopriot in 
LINE 25S 
the lengitudinal mode 
LINE S Se 
Eels 73- page Sa5 


FLOT HEADING 


a0 


LUNE, ts AIRCRAFT WITH AN AUTOPILOT IN 
GBENBS22 =: THE LONGITUDINAL MODE 
EENE 3°: REPe Poe AGE 315 
BOSYOU DESERE TO MAKE ANY CHANGES? 
Ter YES CRY Ne: 
nN 
PUSENG A PRE-ALLOCATED DATASET FOR UNIT -FTI7 F001. 
PUSENG ACPRE-ALLOCATED DATASET FOR-UNIT FII8FoO01. 


A DESS PIA, METARILE HAS BEEN CREATED 
CLEAR SCREEN TO CONTINUE 


ROOT=LOCUS 
OPEN LOOP TRANSFER FUNCTION 
DOL YOU DESIRE TO WAKE ANY"CHANGES- TOs 
1. PLOTPER 7,.PAGE SIZE 
Zee “GRAPH LIMITS 
3. HEADING 
4. GAIN RANGE 
Ss. NO ADDETIONAL CHANGES = PLOT: ROOT=LOCUS 
6. NO-GHANGES — EXIT ROOT=LOCUS PLOTTING 


ROUTINE 
ENTER CPTICN NUMBER. 
6 
RCOT-LOCUS 


OPEN LOOP TRANSFER FUNCTION 
DO YOU DESIRE: 
> “"“GRAPHICAL OUTPUT 
2. TABULAR DATA 
3. OULT SUBROUTINE 
ENTER CPTIICN NUMBER. 
3 
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OPTGR AFH 


OPEN LOOP TRANSFER FUNCTION 


INPUT # = 1 

OQUDEUT 4 =" 4 
DO YCU DESIRE GRAPHICAL RESPONSE AND/OR TABULAR 
CATA FOR: 

Ve. POLE-ZERO MAE 

Z- ROOT=LOCUS 

320 BODE 

4. NYQUIST (POLAR PLOT) 

5. NECHOUS 

6. ANOTHER INPUT/OUTPUT COMBINATION (SAME 

TRANSFER FUNCTION) 
7. ANOTHER TYPE TRANSFER FUNCTION /EXIT 


OPTGRAPH 


ENTER -CPTICN NUMBER. 


5 


OPTGRAFH 


DO YOU DESIRE TO ANALYZE: 


uc 
2e 
Je 
4. 


OPEN LOOP TRANSFER FUNCTION 
NOISE TRANSFER FUNCTION 
COMPENSATOR TRANSFER FUNCTION 
EXIT OPTGRAPH 


ENTER CFIICN NUMBER. 


4 


Re T=12036/ is.tt 20as2e54 
EeCord Off 
END RECORDING OF TERMINAL SESSION 
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The preceding example follows at the end cf the 
GhapGer as figure 3.2 « 


b. Example of tabular data (Terminal Session) 


The entering and exiting option menus for the 
tabular data terminal session are identical to the graphic 
response terminal session option menus, and have been elini- 
nated from the follcwing terminal session. 

BEGENS RECORDING. OF TERMINAL SESSION 
Roe 0-100 102) °20242259 
EXECUTTCN BEGINS... 


ROCT-LOCUS 

OPEN LOOF TRANSFER FUNCTION 
Be YOR DESIRES 

ton” GRAPHTCALOOUIPUT 

2- TABULAR DATA 

Js CULL SUBROUIING 
ENTER CPTICN NUMBER. 
Zz 


RkOOT-LOCUS 
OPEN LOOP TRANSFER FUNCTION 
DO YOU DESIRE) TABULAR OUTPUT WIORGO .TO: 
1.) SCREEN 
Ze. PRINTER 
350, DESKOW(OPGRAPH LISTING) 
ENTER CPTICN NUMBER. 
1 


GAIN KANGE FOR ROOT-LCCUS PLOT 
MINIMUM GAIN = ? 


MAXIMUM GAIN 


N 
oN) 
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GAIN RANGE FOR ROOT-LCCUS PLOT 
MINIMUM GAIN 0.1090D+01 
MAXIMUM GAIN = 0.500D+01 
£0 YOU DESIRE 10 MAKE VANY (CHANGES ¢ 
TYRE UYES® OR: NGM: 


Nn 


HOW MANY PCINTS DO YOU WANT TO TABULATE? 
(500° 1S THE MA REMUM) 
ENTER NUMBER OF POINTS. 


5 
ROOT-LOCUS 
OPEN LOCE TRANSFER FUNCTION 
INPUT NUMBER = 1 
CUTEFUT NUMBER = 1 
TRANSFER FUNCTICN (DC) GAIN = 0.1000D+01 


OPEN LOCF TF POLES ANI ZEROS 


DENCMINATOR ORDER = 
FOLE LOCATIONS 
REAL PART IMAGINARY PART 
0.0000D+00 0.0000D+00 
0.1000D+01 0.0000D+00 
-0.2000D+01 0.3464D+01 
—0.2000D+01 -0.3464D+01 
NUMERATOR ORDER = 1 
ZERO LOCATIONS 
REAL PART TMAGENARY PART 
-0.1000D+01 0.0000D+00 
GAIN = 0.1000L+01 
REAL PART IMAGINARY PART 
0.7074D-01 0.0000D+00 
0.8962D+00 0.0000D+00 
10) a) SSIS) 10) 210) 1 0.3441D+01 
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=O219835D+01 


-0.3441D+01 


GAIN = 0.1800L5+01 
REAL PART IMAGINARY PART 
0.1453D+00 0.0000D+00 
0.7948D+00 0.0000D+00 
=—051970D+01 0.3422D+01 
Oa eOD+:0 4 -0.3422D+01 
GAIN = 0.2600TD+01 
REAL PART IMAGINARY PART 
0.2576D+00 0.0000D+00 
0.6552D+00 0.0000D+00 
=O.1956D404 0.3402D+01 
-0.1956D+01 =-0.3402D+01 
GAIN = 0.3400D+01 
REAL PART IMAGINARY PART 
0.4426D+00 0.1661D+00 
0.4426D+00 =0. 16610400 
-0.1943D+01 O02. 33630D401 
—0.1943D+01 =O. 396s) +,011 
GAIN = 0.4200D+01 
REAL PART IMAGINARY PART 
0.4285D+00 0.3096D+00 
0.4285D+00 -0.3096D+00 
~O-09Z29D+0d 0.3363D+01 
=a 29D +01 =), 33650201 
BO YOU DESIRE TO MAKE ANY CHANGES ? 
TYPESULES VpOR ec" NON, 
ROOT-LOCUS 


OPEN ICCE TRANSFER FUNCTION 
DO YOU DESIRE TO MAKE ANY CHANGES’ TO; 
i. CUTPUT DEVICE 
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2. GAIN RANGE 
3. NUMBER OF POSNTS CALCULATED 
4. NO ADDITIONAL CHANGES - TABULATE DATA 
5. EXIT ROOT-LCCUS TABULAR DATA ROUTINE 
ENTER CFIICN NUMBER. 
5 


ROOT=LOCUS 

OPEN LOOP TRANSFER FUNCTION 
DO YOU DESTRE: 

1. GRAPHICAL OUTPUT 

2- TABULAR DATA 

3. QUIT SUBROUTINE 
ENTER CPTICN NUMBER. 
3 


Re T=T263/3.40 ZOeegesy 
record Gort 
END RECOFDING OF TERMINAL SESSION 


G. FREQUENCY RESPONSE ANALYSIS 


The calculations for the frequency response of the Bode, 
Nyguist, and Nichols responses are identical. The only 
difference is in the graphical presentation of the calcu- 
lated data. The subroutines BODE, NYQST, and NICHOLS cali 
the subroutine FREQ to calculate the transfer function 
Magnitude and phase fcr the frequency range specified by the 
user. 

The subroutine FFEQ eliminates the extraneous zeros by 
comparing the zeros to the Markov parameter as described in 
the Characteristic Equation Analysis section. The frequency 
response magnitude and phase are calculated using double 


precision complex library functions following the 
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conventions of the IBM System /360/370/ FORTRAN IV language. 
The frequency response is calculated by determining the 
product of successive numerator zeros which have been 
combined with the frequency and divided by the preduct of 
successive denominatcr poles which also have been combined 
with the frequency. MThe magnitude for the complex result of 
the zeros, poles, and frequency calculations is determined 
Hsing the FORTRAN |dcubletprecision |) bibrary functions | “for 
finding the magnitude of a complex argument. The phase of 
the of the complex result is determined by separating the 
real and imaginary parts of the complex results, and finding 
the phase angle using the double precision FORTRAN library 
function arctangent. The angle is converted from radians to 
degrees and is corrected for proper quadrant by checking the 
Signs of the real and imaginary parts. The normal range for 
calculating phase angles by the subroutine is from +540 
degrees to -540 degrees. This range is considered adequate 
for most application, but should these limits be exceeded 
Biaes plot) wile sjump with >a’ overtical-linee-to: erther +360 
degrees or -360 degrees respectively before continuing with 
the plot. 

FOE da graphic response the user specified frequency 
Ban ge™ is 9edivided "into =°500 "equally spaced “points (ona 
logrithmic scale) for which to calculate the magnitude and 
phase. The user is given the option to select between 1 and 
500 points for tabular data. 


1. Eode Response 


The graphic response section of the subroutine BODE 
creates two plots. MThe first plot is the response magnitude 
in decibels versus the frequency in ctadians per second. The 
Magnitude from the subroutine FREQ is converted to decibels 
and Plotted with its corresponding frequency. The second 


plot is of the phase in degrees versus frequency in radians 
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per second. The phase from the subroutine FREQ is plotted 
with its corresponding frequency. The frequency for hoth 
plots is plotted .on sa fogrithmicecale.)sThevsea lesrangepror 
both is set automatically to include the entire range of 
data to ke plotted by the subroutine MINMAX which scans the 
data for the maximum and minimum values. System information 
(tranfer function input number, Output number, and (DC) 
gain) is listed in the graph legend. 

The magnitude in decibels and the phase in degrees 
is searched for "0" decibel and 180 degree crossover feints 
to calculate the phase and gain margins. 

The tabular response section gives the user the 
opticn to select the the number of frequency points’ to be 
calculated between 71 and 500 for the frequency range he has 
specified. For takular data response the user is given 
system information (numerator order, denominator order, and 
transfer function (DC) gain) for the selected transfer func- 
tion input number and output number, and the frequency 
(cadians per second), magnitude (decibels), phase (degrees), 


and the real and imaginary frequency response parts. 
2. Bode Response (Interactive Example) 


The fcllcwing example of compensator transfer function with 
a filter and regulator was obtained from [Ref. 8,pp. 382 - 
384}. 

The compensated system is represented by: 


fey = { -254st —19:57 9 
CDT = ee 
[Py = | eon 

| Os =ie6 4 


[On7e7" 1 
GAM} =) olecOare 
15 Oeth 
Ca HP ilosl<oT Ose a) 
[R= b 95.e a 
1 4561. | 

(On Liles 


Example cf Graphic Response (Terminal Session) 


BEGIN RECCRDING OF TERMINAL SESSION 
Reet =0.0:170.02" 18230223 

graphics 

EXECUTION BEGINS... 


OPTGRAPH 


DO TYOU DESIRE TO, AWALYZE's 
se OPEN LOOPH TRANSFER FUNCTION 
Zo RODS: LRANS PERS UNC T LON 
3. COMPENSATOR TRANSFER FUNCTION 
4. EXIT OPTGRAPH 

ENTER OPTION NUMBER. 

3 


OPTGRAPH 
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THE COMPENSATOR TRANSFER SYSTEM OBTALNED FROM 
OPTSYS. CCNDALNS : 


1.. _ NUABER OF STALE. 2 
2 NUMBER OF CONTROLS (INPUTS) = 1 
3. NUMBER OF MEASUREMENTS (OUTPUTS) = 01 


4. MARKOV PARAMETER = 10.0**- 6 


CLEAR SCREEN TO CONTINUE 


OPTGRAPH 


AN OPEN LCOP POLE-ZERO, ROOT LOCUS, BODE, NYQUIST, 
AND/OR NICHOLS PLOT IS DESIRED FOR: 
DNPUT + 2 


CUTEUT “= 2 


OPTGRAPH 


COMPENSATCR TRANSFER FUNCTION 
INPUT # = 1 
OUTPUT # = 17 
DO YOU DESIRE GRAPHICAL RESPONSE AND/OR TABULAR 
DATA FOR: 
1. EOLE-ZERO MAP 
2. ROOT-LOCUS 
3.0 SEODE 
4. NYQUIST (POLAR PLOT) 
5. NICHOLS 
6. ANOTHER INPUT/OUTPUT COMBINATION (SAME 
TRANSFER FUNCTION) 
7. ANOTHER TYPE TRANSFER FUNCTION /EXIT 
CPTGRAPH 
ENTER OPTION NUMBER. 
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BODE PLOT 
COMPENSATCR TRANSFER FUNCTION 
DO YOU DESTRE: 

V1. GRAPHICAL OUTPUT 
2- TABULAR DATA 
, de A UCULT ISUBROUTINE 
ENTER OPTICN NUMBER. 
1 


BODE) PLOT 
COMPENSATCR TRANSFER FUNCTION 
PLOTTER SELECTION AND PAGE SIZE 
(NOTE: PAGE £5 1120 TNCHES BY 8.5 EZNCHES WITH 
SE LECTICN OF TEK6 13) 
WHECH PEOLTER DO YOUCDESERE: 
To TERG 138 
Za VaR oh LEe 
CHOCGSE OPTION 1 OR ©2 


2 
PAGE SIZE (MAXIMUM = 21.0 INCHES BY 21.0 
INCHES) 
HEIGHT = 
5.65 
WIDTH = 
8.5 


PLOTTING (LINTTS TEOR GRAPH 
FREQUENCY RANGE - RADIANS PER SECOND 
MINIMUM FREQUENCY = ? 
0.1 
MAXIMUM FREQUENCY = ? 
1000000 


PEOMTLNG LIMES) IFOR "GRAPH 
FREQUENCY RANGE —- RADIANS PER SECOND 
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MINIMUM FREQUENCY 
MAXIMUM FREQUENCY 


0. 
Or 


DO YOU DESIRE, TO MAKE ANY (CHANGES: 


WE yes 


Nn 


PLOT HEADING 


OR 


wyo" A 


100E+00 
100E+07 


HOW MANY LINESUOF HEADING DO YOU DESIRE ? 


( 3 MAX ) 


ENTER NUMBER OF LINES. 
{ ENTER “O" FOR NO HERDING: ) 


2 
A MAXIMUM OF 32 CHARACTERS PER LINE 15 
ALLCWED 
LINE 1 iS: 
filter simulation 
LINE 2 &S5 
ref. 83 pages: 332 = 334 
PLOT HEADING 
LINE 3 FILTER SIMULATION 
LENE 2 REF 0S? PAGES 332 = 3348 
DO YOU DESIREMTOMMAKE ANY CHANGES 7 
TYPE “YES” OF “HOM. 


A DISSPLA METAFPILE HAS BEEN CREATED 
CLEAR SCREEN TO CONTINUE 


GAIN MARGIN AND PHASE MARGIN 
PHASE CROSSOVER) FREQUENCY iS 75. 188 
GAIN MARGIN IS 
GAIN CROSSOVER, FREQUENCY: 2S) St t43 455 S* 
FHASE MARGIN IS 
CLEAR SCREEN TO CONTINUE 
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=59..8906 “DB 


-89'.9480 DEG 


BODE PLOT 
COMPENSATCR TRANSFER FUNCTION 
DO YOU DESERE TO MAKE ANY CHANGES TO: 
Ht NVELORDER. 7 PAGE SSEZE 
2. GRAPH LIMITS (FREQUENCY RANGE) 
3G “HEADING 
4. NO ADDITIONAL CHANGES - PLOT BODE PLOTS 
5s "NO CHANGES’ = EXIT BODE SUBROUTINE 
ENTER OPTION NUMBER. 
5 


BODE PLOT 
COMPENSATICR TRANSFER FUNCTION 
DO YOU DESIRE: 

Ws GRAPHICAL OUTPUT 
2- TABULAR DATA 
3. QUELT SUBROUTINE 
ENTER OPTICN NUMBER. 
3 


OPTGRAPH 


COMPENSATCR TRANSFER FUNCTION 

INPUT # = 1 

OUTPUT # 1 
DO YOU DESIRE GRAPHICAL RESPONSE AND/OR TABULAR 
DATA FOR: 

1. FOLE-ZERO MAP 

2s) SROOT-LOECUS 

a EODE 

“or  NYQUIST * (POLAR PLOT) 

56 NICHOLS 

6. ANOTHER INPUT/OUTPUT COMBENATION (SAME 
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TRANSFER FUNCTION) 
7. ANOTHER, TYPE TRANSFER: FUNCTION Yer 


CPTGRAPH 
ENTER OPTION NUMBER. 
5 
OPTGRAPH 


DO YOUs DESZRE :TOsANALYZE< 
1. OPEN LOOP TRANSFER FUNCTION 
Ze NOISE TRANSFER) FUNCTION 
3. COMPENSATOR TRANSFER FUNCTION 
4. EXIT OPTGRAPH 

ENTER OPTION NUMBER. 

4 

RR: T=7.0176.604 T82s9219 

record off 

END RECORDING OF TERMINAL SESSION 


The preceding example follows at the end of the 


chapter as: figure 3.3 and “Ezgurce 3.4". 
bh. Example of Tabular Data; Listing File 


The terminal session for the Bode tabular data response is 
Similar to the Bode graphic response terminal session. 
Appendix D contains an example of tabular output sent toa 
disk (OPGRAPH LISTING) for five points between 10 and 100 
radians/second. 


3. Nyguist Respcnse 


The graphic response section of the subroutine NYCST 
creates a polar plot from the magnitude and phase data 
received from the sukroutine FREQ. The plot is automati- 


cally scaled from a minimum radius of one by using the 


51 


subroutine MINMAX to determine the maximum phase and magni- 
tude values. The gain / phase margin calculations and the 


tabular data output are identical to the subroutine BODE. 


4. Nyguist Respcnse (Interactive Example) 


The following example is of the open loop transfer 
function for the filter simulation example described in the 
Bode interactive section. 

This example follows at the end of the chapter as 
figure 325 . 


The graphic response section of the subroutine 
NICHCL creates a single rectangular plot of magnitude in 
decibels versus phase in degrees. The magnitude from the 
subroutine FREQ is converted to decibels and plotted with 
the corresponding magnitude. the plot is) jawtomatically 
scaled using the subroutine MINMAX to determine the maximum 
phase and magnitude values. The gain / phase margin calcu- 
lations and the tabular data output are identical to the 
subroutine BODE. 


6. Nichols Response (Interactive Example) 


The following example is of the open loop transfer 
function for the filter simulation example described in the 
Bode interactive section. 

This example follows at the end of the chapter as 


mgure 3.6 « 
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IV. CONCLUSIONS AND RECOMMENDATIONS 


SS SS Se SS ee Se 


A. CCNCIUSIONS 


An evaluation of the computational ability = GEZOPISYSz 
and CPTIGRAPH, the frogram was tested using an 82 X 82 
matrix, provided by NASA-Edwards, of longitudinal eguations 
for the K-29 experimental forward-swept wing fighter 
aircraft prototype. The OPISYSK program array had to ybe 
redimensioned anda 2-megabyte virtual machine size was 
requized £0r a2 system of this Size. 

The open loop eigenvalues calculated by the OPTSYSX 
program compared favorably with the Eigensysten data 
supplied by NASA-Edwards, but the OPTSYSX program failed in 
calculating the open loop transfer function Zeros. 
Extraneous zeros (calculated by OPTSYSX) were of the same 
approximate magnitude as the open loop system zeros. T has 
prevented the identification and elimination of the extra- 
neous zeros by comparison with the Markov parameter. The 
open loop zero identification was further complicated by not 
having available open loop transfer function zero informa- 
tion from NASA-Edwards to compare the OPTSYSX results to. 

When requested for additional information NASA-Edwards 
Frovided a computer tape with a revised 98 X 98 matrix of 
longitudinal motion equations for the X-29A fighter 
avEcratrt. For this revised system, the OPTSYSX program 
requires a virtual machine size in excess of 2-megabyte. 
Time constraints have prevented running the OPTSYSX program 
with the revised systen. 

The OPTGRAPH program was tested using the OPTSYSX 
program's open loop Eigenvalue and transfer function cutput 
data for the X-29A .82 K 82) j;system of longitudinal motion 
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equations. OPTGRAPH provided excellent quality Bode flcets 
FOr salt of the open loop transfer function inputs and 
eutputs and did not cxhibit any) computational difficulties 
as indicated by the lack of error messages at the end of the 
test. While the test results are encouraging as to the 
ability of the OPTGRAFH program to assist in the analysis of 
large order systems, the actual computational accuracy of 
the OPTGRAPH program can not be fully verified until accu- 
rate transfer functicn information can be obtained from the 
GPToYsSk program. 


Be. RECCMMENDATIONS 


Based on the results of this Thesis the following areas 


emerged as possible areas for further research and study. 
1. Root Finding Routines 


The OPTSYSX program has two similarmroot <iundang 
SubEGutines (HOR and HQR2), to find the transfer function 
Eigensystem (poles) and zeros. The order, of the transfer 
function numerator is calculated separately in the subrou- 
tine ZEROS. : 

The inclusion of the open loop transfer function 
eigensystem (pole locations), zero locations, transfer func- 
tion numerator orders, and transfer function gains, for the 
revised (98 X 98) A1cngitudinal mode system by NASA-Edwards 
should simplify the lcecating the problem in calculating the 
transfer function zeros. 

The open loop system data from NASA-Edwards also 
included Bode magnitude and phase plots, and gain and phase 
Margin information, which can be used to evaluate more 
conclusively the OPTGRAPH program's ability to help evaluate 


a large order system. 
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2. Program Memory Requirements 


The OPTSYSX frogram reguires in excess of 2-megabyte 
of virtual machine memory when configured for large matrix 
operations (98 xX 98). Virtual machines with this remory 
capacity are not normally available toa user. The memory 
usage for matrix storage and manipulation is a possible area 
for the reduction in program memory requirement size. The 
OPTSYSX program retains most of its computational arrays in 
memory. This method simplifies programming but is extremely 


inefficient in its use of virtual machine memory. 


3. Further Modi fications 


The area of modern digital controls should be further inves- 
tigated. The !*computational abilvties of s the OPTCEAPG 
program make it readily adaptable, with minor modifications, 
for assisting in the analysis of transfer functions in the 
digital domains. 
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APPENDIX A 
OPTSYS EXEC LISTING 


EBERACE OFF 
FR RRR RR RK do kkk kkk kk doko tok ak 
THE OPiS ax Ee 
CONTROLS IHE OPTSYSX, OPTCALC AND OPTPLOT 
TOS DETERMINE THE TIME RESPONSE OF A SYSTEM. 
BYeH. XN. DEEL 
VERSION 1.0 16 JULY 1984 


BCDDELED SYM." 1.) LAPT AS 
He KK RK KK RK RK KR KR KK RR RK KR RR RK OK BR OK KOK 


CLRSCEN 
EBEGTYPE -ENDFOUR 


*+ + + + + 


MK KK KKK KKK KK KK He a eK a He He ak a eee eH eK ee a ee ee KK eo ke 


COPPULSFROM OP TSYSK WELLE COME TO THE SCREEN IF YOU 
Wish LaSKh FILE 
ENTER OD ANY OTHER INPUT YIELD SCREEN 


He HH oe He ee He ee eK He a HK ee HK eK KK KK 


-ENDFOUR 

EREAD VAES &ANS 

GIF .&ANS NE .D &GOTO -ONE 

ELLEDEE 06, DISK /OUTPULX, LISTING: Al 
CLRSCRN 

EBEGTYPE -ENDSIX 
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He He ee eK KK ee ee eR eee ee ee eK KK KKK KK KK HK KH 


OUTPUT WILL, GO-TO DISK FILE “OUDPUTK sLesTinG 1A 4 


ee Hee ee Fee ee ee ee it ee eo ok ae ee eae eae eee ae ee ee a ee eo ie ee eK kk ok 


-ENDSIX 

CP SLEEP 3 SEC 

-ONE 

CLRSCRN 

CP LINK 0039P 191 251 RR 

ACGNS1 F 

Ee He ee oe ee ee te he ae eee ae eH ee ee ee ee Ke ea eK eee eK aK KK KK KK 
* * 
* ALLOW THE USE CF AN OLD “OPTMAT DATA Ai" * 
* % 


HR RK RR HR RR RR RR FOR ROR RK Hk OK ROK 
RENAME OPTMAT DATA AT OPTSYS DATA 21 

SIF ERC NEO. €GOTO —FIRst 

RENAME OFTSYS DATA Al OPTMAT DATA At 

FINDSTAK OPTMAT DATA A1 CO1 0 LIMOO02 ALL GROUP1 

ESREAD VARS &A1 &A2 

SIF .&A2 NE .1 &GOTO -START 

CLRSCEN 

EBEGTYPE —ENDTWO 


eR He eR Ha He He a ie ee ee Ree fe he fe ae ae ee ee a ea a eee fe kee he ak a Hee eo oe ok ie ok a oak 


YOU HAVE A DATA FILE NAMED ‘OPTMAT DATA‘ ON YOUR A 
DISK THAT WAS PREVIOUSLY GENERATED BY THE OPTSYsS 
FROGRAM AND CCNTAINS THE F, G, H, GAMMA, A AND B 
MATRICES FROM THAT RUN. 
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rE YOU WOULD LIKE £O USE THESE SAME MATRICES FOR 
TES RUN, THE) €PTSYS PROGRAM WILL READ EN THE 
DESTRED DATA AT THE APPROPRIATE TIME, 


PrP LOU LYPER(Y) ES: 


AN YSOTHER ©ENPUTTWLEL RESULT UN THAT FILE BEEING SERASED! 


FEE EE a Sk OK a He eK ee ee he ee ok ok ek ie ke eo eK ee a 2c ie oe ok a ee KOK OK KK KK KK KK 
-ENDTWO . 
SREAD VARS &ANS 

Ser J6ANS EO, $Y SGOTO =START 


SErRST 
HEE BE HK A ee eK ae He a ae ee 2 ee ek Ke 2K a Fe ie ae He 2 ie ok ee eC aK a Ce a eo ke eo aK Kok 


FRASE THE OLD “OPTMAT DATA Al DATA. FILE 
EEXCE. SOOOP CM SIN THE NEW "OP TMAT DATAOEILE" 
2O ACTIAS APEENG FOR TOPISYSX LAND OPTCALC 


+ + + + + 
FONE > ee eS 


He Re ea aR a RR RK RK KK RK RK RK EK 
ERASE OPIMAT DATA At 

SSTACK 000 0 

FILESTICK OPIMAT DATA Al F 80 1 

=START 

CLRSCEN 

EBEGTYPE -ENDONE 


THE OPTSYS EXEC CONTROLS “FOUR, PROGRAMS): 


© OPTSYSX, FCRIRAN (SYSTEM ANALYSIS) 
GENERATES OPTMAT AND OPTGRAPH DATA SETS 
= RERUN USES) OPTMAT 

Pe IOPTECALGSECRIRAN ~ (CALCULATE? TINECRESPONSE) 
GENERATES OPTPLOT DATA SET USES OPIMAT FROM 
OPTS Vs. 

3 OPTGRAPH FORTRAN (POLE=ZERG,, ROOT-LOCUS, BODE, 
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NYOULS 2. fac) 
USES OPIGROL -OPTGRNO =OPTCRCM FROM OPTot sn 
4 EXIT 


OPTPLOT FORTRAN ES THE FOURTH PROGRAE 
USES OFIPLOT DATA SET FROM OPTCALEC 
TEE SIZE OF. THE DATA SETS: VARA -WLIE SIRE Sia 2eu 
ORDER, AND CAN USE ABOUT 20% OF THE USERS SOLS hese ser. 
THEREFORE ENSURE THAT SUPFFICLENT DISK SExGE ts 


AVAILABLE. 
IF DATA SET IS ALREADY AVAILABLE YOU MAY RUN ANY 
PROGRAM 
= ENTERS Si geen gees eet ANY OTHER ENEUT 
RETURNS TO MENU 
—-EN DONE 


SREAD VAES 6&ANS 

SLF .SANS EQ «1 GGOTO —OPTS¥sS 

SIF .€RNS EQ 22 EGOTO —OPTCALC! 

SIF .&ANS EQ .3 &GOTO -—OPGRAPH! 

Si PF GANS? EO «4 GCEXLT CRE 

&GOTO -START 

=OPTSYs = 

FILEDEF <8 DESK OPTPLCT DATA AT s(BERM 
FIGEDER "9 DESK (OPEMAT DATA AW (CER 
FILEDEF 0 DESK OPTEGROL DATA Al (PERM 
FILEDEF 01 DISK OPTGRNO: DATA Al (PERM 
FILEDEF O04 DISK OPTGRCM DATA Al (PERM 
GLOBAL TXTLIB VFORTLIE IMSLDP NONIMSL 
ETYPE LOADING OPTSYS.... GENERAL STATE VARIABLE ANALYSIS 


PROGRAM 
GLY PE craton OPTMAT DATA Al w.s0- INPUT DATA SET CNP REAUNS 
OPTS 1 SK 
CLRESCEN 
STYP Ere) DADA SES” (Or sana n AND .- OPGRAPH ‘CREATED 


CP SLEEP 5 SEC 
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GOTO" sSSTART 


SGPTCALTE1 
SE Se RK ae RR a KK KKK KKK KK RK OK RK 


* x 
= CHECK FOR DATA IN THE FILE "OPTMAT DATA " * 
* BEFORE LOADING OPTCALC * 
* * 


BRE He EE ee ee He 2 ee ee ee ie 2 2c ie ie ok oe ok ok ae oe 2 ok ok ok OK ie Ko KK KK KK KK 
RENDSTARK LOPTIMAT DATA ‘A1°-CO1 O° LIMOO2 ALL GROUP 
SREAD VARS &Al1 &A2 

SEE each GEO's EGOTO -=OPTCALC 

SrYPETPROPER DATAVFILE 15 NOT AVALDABLE FOR. OPTCALC 
EPPSLEED(S SEC 

&GOTO -START 

=GRECAEC 

GEYPE.ssse0e56s0PTCALC.... \TEME RESPONSE PROGRAM 
ELGEDER DE (DUSK OPTPIUCT DATA Al (PERM 

PEEEDEF 9 DISK OPTMAT DATA Al (PERN 

GLOBAL TXTLIB VFORTLIE IMSLDP NONIMSL 


OPTCAIC 

STYPE oO PMPLOT, DATA tA CREATED 
GEBTSEEEP 5 SEC 

CLRSCRN 


EBEGTYPE —-ENDNINE 


eK ee eo ee He eK ee ee ee eK a Ke He ee ae ee ee OK IK 


TOU ARE ehISsSh LEST PED Wil THE RESULTS 
PAUSY FAR SAND WOULD ULEKESTOMExXET ‘ToOr'Cns, 


VP toe Or ber 
(ANY POTHER NPUTYTONCONEENUE) 
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MEK KK a eK I eK a a oe ee eee He ee ie fe ee a 2 ae ee ke a a ee i ke ek aK ok ok 


-ENDNINE 

GSREAD VARS 6ANS 

GLF .GANS EQ .ieCOTO —STaAra 

GTYPE FLOTS NEXT 

GP SLEEP &@ SEC 

GLOBAL TXTLIB DISLIBVS 92DISETBR S2ZINTLIB VFORTLIBD ”*GRFELTE 
NONIMSL CMSLIB 

FELEDEF 6 DISK OPTPLOYT DATA Al (PERM 

FILEDEF 11 DISK LBLANK TMP A3 (RECREM PF TREC 2400 BL eSiLZe 
2400 BLKSIZE 2400 XTENT 600 

FILEDEF 12 DISK ISSCCMAP MAPDTA * (RECPM F TREEr “0G 
BLKSIZE 400 

FILEDEF 13 DISK TABLET TMP A3 (RECFM VS LRECL 208 BLURS 
208 XTENT 1000 

FILEDEF 14 DISK LBLANK MAPDTA &LBMODE (RECFM VS LRECL 608 
BLKSIZE 608 

FILEDEF 17 DISK DISSFPLA SYSsuT1 AS e(RECEM EER ECL 2000 
BLCCK 2000 XTENT 500 

FILEDEF 18 DISK DISSELAVMETAFILE AS. ((RECEMOOVESELEECL 191Ges 
BLCCK 19069 


OPTPLOT 

&GOTO -META 

-OPGRAPH1 

Ae fe ee ae eK eae ee eae ae ae ie ae Se eo a he ae fe he ee ee oe a ee ee KK ie ee He a Ke eK KK Kk 
xk * 
* CHECK FOR DATA IN THE OPTGRAPH DATA SETS " = 
* BEFORE LOADING OPGRAPH * 
x = 


HR HC HR RK HR HK KK KR RR OR RR I KR OK ak 
SEH SE Fe HR Ro I RR SR eo eR RR OK RK OK RO OK EK 
RENAME OFTGROL DATA A1 OPTSYS DATA Al 

SIF ERC NE 0 EGOTO =CETNOISE 

RENAME OPTSYS DATA A1 OPTGROL DATA Al 
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FINDSTAK OPTGROL DATA A1 CO1 0 LIM002 ALL GROUP1 

SREAD VARS &A1 &A2 &SA3 GA4 SAS 

GIF .6A2 NE .0 §GOTO -SIXTEEN 

-OPINCISE 

RENAME OFTGRNO DATA A1 OPTSYS DATA Al 

SIF &RC NE 0 &GOTO -CETCOMP 

RENAME OPTSYS DATA A1 OPTGRNO DATA Al 

FINDSTAK OPTGRNO DATA A1 C01 0 LIMOO2 ALL GROUP 

G6READ VARS 6A1 &A2 &A3 5A4 €A5 

SIF .&A2 NE .0 &GOTO -SIXTEEN 

-OPTCCMP 

RENAME OPTGRCM DATA A1 OPTSYS DATA Al 

GIF &RC KE 0 &GOTO -NCDATA 

RENAME OPTSYS DATA A1 OPTGRCM DATA A1 

FINDSTAK OPTGRCM DATA A1 CO1 0 LIM002 ALL GROUP1 

SREAD VARS &A1 &A2 &A3 SA4 EAS 

SIF .6A2 NE .0 &GOTO -SIXTEEN 

-NODATA 

SPRINT PROPER DATA FILE IS NOT AVAILABLE FOR OPGRAPH 

CP SLEEP 5 SEC 

&GOTO -START 

-SIXTEEN 

STYPE ...OPTGRAPH.... CLASSICAL ANALYSIS OF OPTSYS OUTPUT 

SET LDRTELS 10 

GLOBAL TXTLIB DISLIBVS 92DISLIB 92INTLIB VFORTLIB GRFLIB 
NONIMSL CMSLIB 

SET LDRTBLS 10 

FILEDEF 01 DISK OPTGRNO DATA A1 (PERM 

FILEDEF 02 DISK OPGRAFH LISTING A1 

FILEDEF 03 PRINTER (RECFM FA BLKSIZE 133 PERM 

FILEDEF 04 DISK OPTGRCM DATA Al (PERM 

FILEDEF 06 TERM (RECFM FA BLKSIZE 133 

FILEDEF 05 TERM (RECFM FA BLKSIZE 80 

FILEDEF 10 DISK OPTGECL DATA A1 (PERM 

FILEDEF 11 DISK LBLANK TMP A3 (RECFM F LRECL 2400 BLKSIZE 
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2400 XTENT 600 

FILEDEF 12° DISK ISSCCMAP MAPDTAV> CREECH FeOR EGE 400 
BLKSIZE 400 

FILEDEF. 13 DISK TABLEI TMP A3 (RECFM VS LRECL 206 SER auq 
208 XTENT 1000 

FILEDEF 14 DISK LBLANK MAPDTA ELBMODE (RECFM VS LRECL G08 
BLKSIZE 608 

FILEDEF’ 17 DISK DISSELA SYSUT? AS (RECEM FR eRRECT (2000 
BLOCK 2000 XTENT 500 

FILEDEF 18° DISK DISSPLA: METAFI LErAarr(RECPNNVES) LRECE si 0G4 
BLECK 19069 


OPTGRAPH 

-META 

toe qoackgackduck 4 dick 44 idekemerkerd anekdcaaane 840k ears a 
3 *% 
* CHECK FOR FILE "DISSPLA METAFILE A4" ON * 
% THE USER'S DISK BEFORE GOING TO DISSPOP = 
* * 


He KK RK Ro oo oo oo a a aK oe oko fk ee ok KK kK ek 
RENAME DISSPLA) METAFPALIE AY OPTSYS METAFILE AG 

SIF GRETECODE NE OVEGETO =STARZ 

RENAME OPTSYS METAFILE A4Y DISSPLA METAFILE A4 

=ELGHT 

CLRSCEN 

SBEGTYPE —ENDTEN 


HR I He ee a 2 He ee ok Hee ee oe He ke ee 2 a HK eK eK ke Ke KO He KK kK KK 


DC: YOU WANT A VRSZEC PLORTER] SMOOTH COPY OR TTHE 
THE. DISSPLA, METAETLE THAT YOU RIUSH, CREATED? 
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(Y OR N) 


FR eR 2 ee ae a ae ee a a ae ke eo ko oe oo oo ok eo ok ok oko ok ok oe ok Ko ok kek oe ek kk 


-ENDTEN 

SREAD VARS &ANS 

GEF =~cGANS EQ .Y¥ “SGOTC “NINTH 
cere] CANS: EQ) oNseGOTC —START 
&GOTO -EIGHT 

-NINTH 

FEEC.DISSPOP VRSTEC 

CLRSCEN 

SBEGTYPE -ENDTWELVE 


FETE EK OK A HK HK eK eK ee a ak eK a EE aK ee Ko ok ae kk 


FOUR CRAPH(S)- CAN -BE| PICKED UP AT -THE. COMPUTER -CENTER. 
HE. GRAPH (SS), WLLL BE. ADDRESSED, -1TO_ "POR: (USER LD)". 


WRK HE eK eK KK He a ee ee ee He 2 Se a 2 He ee A eK I KO EK ok KKK KK ee 


-ENDTWELVE 
CE SLEEP 5 SEC 
&GOTO -START 
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APPENDIX B 
OPTSYSX PROGRAM LISTING 


This portion of the thesis contains the OPTSYSX FORTRAN 
program (88 pages). 
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